Polarization Measurements and the Pairing Gap in the Universal Regime 



J. Carlson and Sanjay Reddy 
Theoretical Division, Los Alamos National Laboratory Los Alamos, NM 87545, USA 

(Dated: February 2, 2008) 

We analyze recent cold-atom experiments on imbalanced Fermi systems using a minimal model in 
which a BCS-like superfluid phase coexists with a normal phase. This model is used to extract the 
T=0 pairing gap in the fully paired superfluid state. The recently measured particle-density profiles 
are in good agreement with the theoretical predictions obtained for the universal parameters from 
previous Quantum Monte Carlo (QMC) calculations. We find that the zero-temperature pairing gap 
is greater than 0.4 times the Fermi energy Ef, with a preferred value of 0.45 ± 0.05 Ef- The ratio 
of the pairing gap to the Fermi Energy is larger here than in any other system of strongly-paired 
fermions in which individual pairs are unbound. 
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Experiments that trap and cool fermionic atoms are 
now providing new insights on the nature of strongly- 
interacting Fermi systems with a number imbalance be- 
tween the interacting species [J, [3, H, 0, [1| ■ Experiments 
to date have studied systems containing two hyperfine 
states of ^Li, which we label | t) and | |) for convenience. 
These experiments can tune both the number asymmetry 
(polarization) and the interaction strength and thus have 
the potential to probe new phases of superfluid matter 
that are expected on theoretical grounds. Furthermore, 
they can directly probe pairing gaps in these strongly- 
paired superfluid systems. In this paper we examine re- 
cent experimental data on polarized Fermi systems and 
extract the T=0 pairing gap in the BCS phase. The ratio 
6 of the pairing gap to the Fermi energy at unitarity is a 
universal parameter and previous QMC calculations pre- 
dict 5 « 0.5. Gaps in strongly-paired Fermi systems are 
important in diverse areas of physics including nuclear 
matter in neutron stars and the phase structure of dense 
QCD. 

Cold-atom experiments probe the strongly-interacting 
regime where the superfluid properties are robust at fi- 
nite temperature. When the short-range interaction is 
tuned to produce an infinite scattering length all measur- 
able quantities are related to their free Fermi gas coun- 
terparts by universal constants since the interaction does 
not present a dimensionful scale. For this reason, the sys- 
tem with a = cxD is said to be in the "universal" regime. 
In this regime two phases are certain to exist: a super- 
fluid state at zero polarization and a normal state at large 
polarization. The superfluid state at T=0 can be char- 
acterized by the ground-state energy Esf = ^(S/S)^^^ 
and the pairing gap A = SEp with Ep = (37r^p)^/'^/2m. 
The normal state can be characterized by the binding 
of a minority spin particle to the Fermi sea of majority 
particles, En — Eq = —xEf^N, where En is the energy 
of the normal state, Eq is the energy of non-interacting 
particles, and Ef^N is the Fermi energy of the majority 
spin population Ef^N — p)^^^ /2m. 

The universal constants have been calculated using 



Quantum Monte Carlo techniques, @, 0, B B yield- 
ing: ^ = 0.42(01), d = 0.50(03), and x = 0.60(.01). 
More recent calculations of ^ using both Diffusion Monte 
Carlo and Auxiliary Field Monte Carlo [lH indicate that 
^ may be slightly smaller, ^ = 0.40(.01). While these 
calculations are approximate, calculations of the respec- 
tive ground-state energies are upper bounds and hence 
provide an apparently accurate upper bound to the pa- 
rameter f and lower bound to x UM- The pairing gap 
parameter 6 is obtained from the difference between even- 
and odd-particle number simulations, and hence is not a 
bound on the true value. Measurements of the parame- 
ters ^ and X appear to be consistent within errors with 
these calculations as shown below. 

These universal parameters are deflned precisely at po- 
larization zero and one for the superfluid and normal 
states, respectively. QMC calculations of the superfluid 
and normal state suggest that simple descriptions based 
upon a quasiparticle picture and these universal parame- 
ters are valid over a wide range of polarizations between 
zero and one. At T = the superfluid -I- quasiparti- 
cle picture appears to work reliably up to polarization 
a = (n| — ni)/{n-] + nj.) < 0.2 while in the normal 
state the simple quasiparticle picture appears to work re- 
liably for polarizations a > 0.4 [lol]. The phase transition 
between these two states is flrst-order at low tempera- 
ture. At T=0 the transition occurs between an unpolar- 
ized superfluid state and a normal state at a polarization 
of approximately 0.4 [lo| . At the transition — dor u 
where Sfi = (/i| — fii)/2 and fj, ^ {fi^- + fJ.i)/2 [9|, llOj . 
Here the spin-up and spin-down particle chemical poten- 
tials are /if and /i| , respectively. The pairing gap is only 
slightly larger than the value required for the stability of 
a gapless superfluid with finite polarization at zero tem- 
perature. At finite temperature, however, the superfiuid 
quasiparticles will be thermally excited, polarizing the 
superfiuid; this polarization can be exploited to extract 
the pairing gap in the superfiuid. The main goal of this 
study is to extract these fundamental properties of the 
system from experiments in the universal regime. 
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The first-order phase transition appears to have been 
observed in recent MIT Q and Rice Q experiments. Sig- 
nificant differences remain, however. In particular the 
Rice experiment finds a transition directly from zero to 
full polarization. Such a scenario cannot happen in a 
bulk system if the binding of spin-down particles in the 
normal state (x = 0.6) is so large, and hence the observed 
polarization must be related to the unique geometry of 
the trap. The MIT experiment has been performed for 
larger numbers of particles and more spherical traps, and 
hence the simple local density approximation employed 
here is better suited for an analysis of these experiments. 

Our picture of the two states at very low temperatures 
are presented below. The calculated universal parame- 
ters then yield precise density distributions for the mi- 
nority and majority species in the trap. 

Superfluid state with a small polarization: In the uni- 
versal regime the energy density, pressure and the chemi- 
cal potential of the unpolarized superfluid state are given 
by esF = C A:|/(107r2TO), Psf = C kl/{lbTr^m) and 
/i = ^ kp/2m, respectively. Here /cf is the Fermi mo- 
mentum kp — (i-K^pf-l^ . In Ref. Q we calculated the 
excited state quasi-particle (qp) dispersion relation by 
introducing additional fermions to the unpolarized su- 
perfluid state. To leading order in the momentum ex- 
pansion, the dispersion curve for qp's measured relative 
to /i is given by 



a;qp(fc) = Av/(1 + a2 (a: - xo)^ + 0\{x - xo)^)] , (1) 

where x = feV^F, a^o ^ 0.83(3), and a-i ~ 1.3(2). The 
higher order terms are negligible at small polarization. 

The energy to excite a quasiparticle in the vicinity of 
the normal-superfluid transition is anomalously small in 
strong coupling since S^c ^ A at the first-order transi- 
tion [9]. Thus even when T ^ Ep spin- up quasiparticles 
are thermally excited in the superfluid in the vicinity of 
the transition. Further, zero-temperature QMC calcula- 
tions suggests that interactions between quasi-particles 
in the superfluid are weak at low polarization. Hence the 
number density and pressure of thermally excited quasi- 
particles are taken to be 
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respectively. The superfluid phonons with a dispersion 
relation lu — c k, where c = kp/rn, are negligible 

compared to spin-up quasi-particles except at very low 
temperature when T ^ (A — (5/i). 

We anticipate that the gap will decrease due to both 
finite temperature and polarization. For the present 
calculations we assume that for low temperature the 
pairing gap decreases slowly according to A(T) = 
A(0)^1 — {T/tc EpY). The zero-temperature pairing 



gap is a constant fraction 5 of the Fermi energy: A(0) = 
5Ep. The coefficient tc controls the temperature depen- 
dence of the gap and depends on A(0) and the ratio 
Its calculation is beyond the scope of this work and a 
mean-field treatment of the polarized finite temperature 
superfluid phase is discussed in Ref.[l3|]. Here we treat 
tc as a parameter in the range 0.05 — 0.25. 

Normal state at high polarization: Through QMC 
studies we have determined the energetics of the normal 
state at high polarization. Near unit polarization, earlier 
work has shown that a mean field shift for the spin-down 
particle given by 
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with X — 0.6 provides a good description of QMC re- 
sults for the energy at zero temperature [iq| . To ex- 
tend the analysis to finite temperature we adopt a sim- 
ple independent-particle model in which the interactions 
modify the single-particle levels and are able to fit the 
QMC results. A symmetric form for the dispersion rela- 
tion that fits the QMC results is given by 
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where kp = {k%^ + k%^Y/^ and x = 0.6. We note that 
the dispersion relation in Eq. 3] and the symmetric form 
in Eq. [5] both provide a fair description of the QMC data 
near unit polarization and differ only by a few percent at 
cr = 0.5. 

We assume a harmonic trap and work in the local 
density or Thomas-Fermi approximation. The physical 
state at any location in the trap is completely deter- 
mined by the local chemical potentials /i| = fi + Sfi and 
111 — ^ — 5 PL. The chemical potential /i = A — VTrap('^), 
where VTrap('') = 0.5 hio (r/ro)'^ and ro = hc/Vhwrru? 
while 6)1 is constant since the trap does not distinguish 
between the different fermion species. Using the notation 
of reference distances are scaled relative to the radius 
i?l at which the zero-temperature spin-up density van- 
ishes. The radius is similarly defined as the radius at 
which the zero-temperature spin-down density vanishes. 
The temperature T' is scaled to the Fermi energy of a 
non-interacting Fermi gas at r = with the same den- 
sity profile as the exterior cloud. Similarly the densities 
are measured in units of no, where no is the density of 
a non-interacting Fermi gas at r = with the same den- 
sity profile as the exterior cloud. At low temperature the 
density profiles predicted by the theory discussed above 
are shown in Fig. [TJ The inner superfluid region and the 
outer normal region separated by a first-order transition 
are clearly seen. The transition is characterized by as 
and cr„ - the polarizations in the superfluid and normal 
phase, respectively. As mentioned earlier, at low tem- 
perature ajsr ~ 0.4 while (7^ increases rapidly from zero 
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FIG. 1: Spin-up and spin-down densities and polarization 
versus radius predicted by theory for 5 = 0.43 at T' = 
(dashed) and T' = 0.03 (sohd) are shown. The universal 
parameters ^ and x s^rs taken as 0.4 and 0.6, respectively. 
For T' — 0.03 the local value of S shown by the dot-dashed 
curve was obtained using tc = 0.1. 

at zero-temperature with increasing temperature due to 
the low excitation-energy of spin-up quasiparticles in the 
vicinity of the transition. In our implementation of the 
equation of state for the normal phase, thermal effects on 
the density profiles are negligible except close to r = Ri 
and r = i?| . 

The density distribution of the spin-up particles (ma- 
jority species) between and i?| provides a measure of 
the temperature and the chemical potential of the spin- 
up particles. Similarly, the density distribution of the 
unpolarized superfluid state at the origin provides an in- 
dependent measure of the average chemical potential 
assuming the calculated value of ^. Since x = 0.6 de- 
scribes the normal state near unit polarization, the ra- 
dius Ri can independently provide a measure of the spin- 
down chemical potential. Thus, the two radii, the density 
at large radii, and the central density measure the two 
chemical potentials and the temperature and provide a 
consistency check between the calculated values of ^ and 
X in extracting ^ and Sfi. 

In order to compare with experimental results, we as- 
sume a specific ratio of the density at the origin to the 
maximum density of a cloud of free fermions with the 
same density distribution for r > which is equiva- 
lent to assuming a specific average chemical potential /i. 
The values of the ratio RJR^, the total polarization Ptot 
in the trap, the transition radius i?c, and the polariza- 
tion as a function of radius are then predicted for various 
values of the pairing gap and the temperature. We find 
that present measurements provide strong constraints on 
the pairing gap even though the temperature is not very 
precisely determined in the lowest temperature measure- 



FIG. 2: Spin-up and spin-down densities theory and experi- 
ment for 5 = 0.43, 5 = 0.4 and x = 0.6 at T' = 0.03. 

ments. 

Fig. [2] compares the calculated spin-up and spin-down 
densities as a function of radius. We used a normalized 
temperature T' = 0.03 and a normalized density at the 
origin n-\{r = 0)/no = ni{r = 0)/no = 1.72 consistent 
with the experiment. This simple model reproduces the 
radius of the transition, the radius where the spin-down 
density goes to zero, and the overall polarization in the 
trap. The calculated polarization is 0.44 for these values 
of the parameters, and the measured value is 0.44(0.04) 

In Fig. [3] we compare the measured and calculated 
polarization at T = 0.03 and T = 0.05. Results for 
different values of 5 and tc are shown. At T' = 0.03 
the polarization is approximately 0.12 at the interface 
at r/_R| — 0.43, consistent with the experimental results. 
The qualitative and quantitative features of the measured 
polarization at T' ~ 0.03 are captured by the normal 
phase at r/R^-^ 0.45 and a thermally polarized super- 
fluid phase for r/R^ < 0.4 - consistent with a first-order 
transition somewhere in between. In contrast, the com- 
parison between theory and data at T' — 0.05 suggests 
that the superfluid extends further out. Extrapolation 
of the polarization predicted by theory in the superfluid 
state (dotted curve) to p ~ 0.4 provides a much better 
description of the data than the normal state prediction 
and a clear signature of a first-order transition is absent. 
In both cases there appears to be evidence for a mild de- 
crease in the gap with increasing T/Ep and polarization. 

For a fixed central density and i?| , our minimal model 
predicts that the phase-boundary Rc moves outward in 
the trap with increasing temperature. This behavior is 
very sensitive to the thermal properties of both phases 
at low temperature. In our analysis, which is expected 
to hold only for small temperature and polarization, the 
thermal response of the superfluid phase in the vicinity of 
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FIG. 3: Polarization versus radius, theory and experiment, for 
different values of 5 and tc at T' = 0.03 and T' = 0.05. The 
dashed curves show the local finite temperature gap. The 
results indicate that the data provide both an upper and a 
lower bound on the gap: 0.5 > (5 > 0.4. 



the transition is stronger than that of the normal phase 
- driven entirely by the fact that spin-up quasiparticles 
are easy to excite and have a large density of states. 

The comparison in Fig. [3] provides compelling lower 
and upper bounds for the superfluid gap. Even if the 
temperature was extracted incorrectly from the exper- 
iment, the extracted gap cannot be too small. A gap 
smaller than « ^AEp would produce a shell of polarized 
superfluid before the transition even at zero temperature. 
Furthermore, the radial dependence of this polarization 
would be quite difl'erent than observed experimentally, 
rising abruptly from the point where A = 5/i and being 
concave rather than convex. A gap larger than w 0.5-Ei? 
would be unable to produce the observed polarization in 
the superfluid phase. 

We have also examined the dependence of our results 



on the universal parameters ^ and x- Both of these are 
expected to be uncertain by 0.02. These uncertainties, 
as well as the uncertainties in the superfluid quasiparticle 
dispersion relation do not significantly alter the extracted 
bounds on the superfluid gap. If the transition tempera- 
ture and gap decrease significantly at finite temperature 
with increasing polarization, the extracted gap could be 
slightly higher. 

We note that the pairing gap extracted here is signifi- 
cantly larger than that obtained in simple analysis of RF 
spectroscopy. These experiments measure the response 
to a probe tuned to a transition between a mjnoritj 
for example, and a third hyperfine state 
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rity spin. 



RF spectroscopy analysis, however, is complicated by the 
strong final-state interactions that reduce the average en- 
ergies of the transition[l7l|. 

In summary, it is possible to extract the pairing gap 
from measurements of polarized Fermi gases in the uni- 
tary regime. These systems have an extremely large gap 
of almost one-half the Fermi energy - the value extracted 
in this work is clearly the largest gap measured in any 
Fermi system. Precise measurements of the particle den- 
sities in the unitary regime, and measurements extending 
into the BCS regime would be valuable, in particular as 
direct experimental tests of the pairing gap in low-density 
neutron matter which is relevant in neutron stars. 
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